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Abstract 

The results of Brownian dynamics simulations of a single DNA molecule in shear 
flow are presented taking into account the effect of internal viscosity. The dissipa- 
tive mechanism of internal viscosity is proved necessary in the research of DNA 
dynamics. A stochastic model is derived on the basis of the balance equation for 
forces acting on the chain. The Euler method is applied to the solution of the model. 
The extensions of DNA molecules for different Weissenberg numbers are analyzed. 
Comparison with the experimental results available in the literature is carried out 
to estimate the contribution of the effect of internal viscosity. 
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1 Introduction 



Along with the progress in study of biological properties of polymer (such as 
DNA) molecule s, their mechanical features in solutions have also been inves- 



tigated recently (ISmith et al.l . Il999l ). Modeling the macromolecules is a key in 
studying the dynamics of single molecules or solutions. The bead-spring chain 
model is a widely used model to simulate the chain molecules. If the entropic 
springs connecting the beads are assumed nonlinear, the governing equations 
of the system can not be solved analytically. The development of efficient nu- 
merical procedures becomes then necessary. Indeed, the increasing numbers 
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of beads in the molecule chain makes the numerical simulation difficult and 
time-consuming. However, the dumbbell model with only two beads is known 
to be computationally efficient. Furthermore, it can yield good results if ap- 
propriate factors, such as internal viscosity, hydrodynamic interactions, are 
considered. 



There are two main approaches in investigating pol ymers in solut ions: one 



is to study the bulk rheology of polymeric solutions (ILarson 1. 120051) and th e 



other is to chase the dynamics of a single molecules in detail ( IShaqfeh 1 . 120051 ) 



With the ability to examine dynamics of macromolecule chains by fluores- 
cence microscopy measurement developed in recent years, the research on the 
deformation and dynamics of single molecules have progressed substantially. 



In the study of bulk rheology of polymer solutions, the dynam ics of the flow of 



start - up or sudden cessation h ave been discussed t heoretically (IPrabhakar et al. 



2005 



2002 



Hua fc Schieber 



Goshen et al. 



■..19951) and experimentally (jSmith fc Chul . Il998l ; iGoff et al. 
2005f ). [Hua k, Schieber I (Il995l ) studied the finitely exten- 



sible nonlinear elastic (FENE) dumbbell model with internal viscosity (IV). 
They discussed the contribution of IV to the bulk rheology of polymer solu- 
tions. 



On th e other hand, by video fluorescence microcopy method, iSmith et al. 



( 119991 ) experimentally studied the fractional extension of A-phage DNA molecules 
in steady shear flow for different Weissenberg numbers that denote the dimen- 
sionless shear rates. They also calculated the probability distribution of the 
extension by a statistical method. Although the bulk rheology can be obtained 
by further computation of the extension of the dumbbell model molecules, the 
analysis of the extension and its distribution for s pecifi c polymers analytically 



is still limited in the literature flHur fc Shaqfeh LboOOh . 



In this paper, the FENE dumbbell model with IV proposed by lHua fc Schieber 



( 119951 ) is employed to study the extensions of A-phage DNA molecules in 
steady shear flow. After numerical computation, the results of molecule ex- 
tension a nd probability distr ibution are compared to that obtained in exper- 
iment by ISmith et al.l (Il999l ). The plot of fractional extension for different 
Weissenberg numbers is also superimposed by that of different models, e.g., 
the bead-spring model with wormli ke spring bvlJendreiack et al. I (120021 ) and 



classical FENE dumbbell model by iHur &: Shaqfeh I (120001 ) . We show that the 
FENE dumbbell model with IV can yield good agreement in simulation of 
A-phage DNA molecules. 



In Section 2, we give a brief description of the dumbbell model with IV and 
its governing stochastic differential equation. In Section 3, the details of the 
algorithm to solve the governing equation and the values of parameters are 
given. The numerical results are presented and discussed in Section 4, followed 



2 



by concluding remarks. 



2 Model for the Dynamics of DNA molecules 

From classical physics viewpoint, the mass points, beads of the dumbbell 
model are assumed to be subjected to three forces: 1) a hydrodynamic force, 
which describes the resistance to motion of the bead bathing in the solvent, 
2) a connector force, which is denoted by the spring between the beads, 3) a 
Brownian force, which is caused by the rapidly bombarding of the beads by 
small solvent molecules. The second force is due to the entropic force leading 
the molecules to their statistical equilibrium positions. It is natural to think 
that there exists, a forth force, the resistance to motion experienced during 
changes in the molecule conformation. In what follows, the internal viscosity 
force is introduced to reflect such intramolecular resistance effect. 



Figure 1. The dumbbell model with IV Figure 2. The shear flow 

We start our consideration with the bead-spring-bead dumbbell model of a 
single DNA molecule in a Newtonian solvent with viscosity i] s . It is assumed 
that there is no interaction among the beads of different dumbbells and the 
viscous drag coefficient is due to the resistance of the flow denoted further by 



For the dumbbell model with internal viscosity, we use a schematic representa- 
tion of Figure 1 where Q is the connector vector of the two beads. The spring 
force is a function of the configuration vector and configuration velocity. This 
force is the combination of connector force and the IV force mentioned in the 
first paragraph of this section. If we consider the finitely extensible model, we 
get a force law of the following form: 



where extension Q defined by ^/Q • Q. In equation (1), Q is the maximum 
spring extension, and when the extension of the dumbbell approaches this 
value, the spring force tends to be infinite; H and K are Hookean constant 




X 




(1) 
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and interval viscous coefficient, respectively. The second moment tensor QQ 
is a key to study the dynamics of the molecules. The first term of the ri g ht sid e 
of equation (1) denotes the FENE spring force proposed by IWarnerl (119721 ). 
Considering the forces described above, the equation of one bead's motion 
can be ob tained if the inert ial term is neglected. According to the phase-space 
theorem (IBird et all 119871 ). the equation of motion can be given in the form 
of the equation of the dumbbell configuration vector Q: 




(C/2JT) + 1 Q 2 




2kT d 2 
k, ■ QJ — -— ^ In if) — — 



H 



C 1 - (Q/Qo) 



rQ 



(2) 

where ib is the configuration distribution function of Q changing with time. 
The fluid velocity is given by specifying transpose of imposed fluid velocity 



gradient k, 



fVv) 



and d is unit matrix. 



Now we introduce dimensionless parameters 



Q 



Q 



AHt 



kT/H 



- t 



(3) 



so that equation (2) could be cast into the dimensionless form 



dQ 

di 



^=5 



QQ 



(C/2K) + 1 



\k ■ Q 



- — =- In ib 

2dQ V 2 



Q 



- \ 2 



Q Qo 



(4) 

where A = (/4:H is the relaxation time of the molecules and \JkT /H is the 
root-mean-square average size of Hookean dumbbell in one dimension at equi- 
librium. For convenience, the bars over the dimensionless parameters are omit- 
ted thereafter without confusions. 



Before we get the governing stochastic differential equation, the Smoluchowski 
equation for this dumbbell model should be deduced, which can be obt ained 
by substituting equation (4) in equation of continuity (IBird et all 1 198 71 ): 



dip 
~dt 







dQ 




(5) 



where the internal viscosity parameter is defined as e = 2K/(. 
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With the Ito interpretation we can obtain the equivalent sto chastic differential 
equation of equation (5) (see, e.g.. lHua fc Schieber l jl995h ): 



dQ 



QQ 



e + 1 Q 2 
+ 



1 Q 


\_ 


e Q 


21-(Q/Q )V 


e + lQ 2 


1 \ QQ 


■dW t , 




s + l) Q 2 





dt 



(6) 



where 



(dW t ) = 0, 

(dW t dW t ) =6(t- f) 5dt. 



(7) 



If e = in equation (6), the stochastic differential equation de generates to case 
for th e clas sical FENE dum bbell case as in the references by iHur fc Shaqfeh 
(l2000h and lOttingerl dl996h . 



In the next section, the algorithm and the specific parameter values will be 
presented in detail for A-phage DNA in steady shear flow. 



3 Computational Implementation 



The steady state simple shear flow in the Cartesian coordinate system has the 
velocity vector v in three dimensions, with elements v x = 'jy, v y = 0, v z = 
and where 7 = dv x /dy is the shear rate (See Figure 2). In such flow, k can be 
expressed in linear function of shear rate 7 



^0 7 0^ 








So we can adjust the parameter \k in equation (6) by tuning Weissenberg 
number Wi defined by W{ = Xj. 

Based on Euler's method, equation (6) can be discretized as: 



Q 



t+At 



Qt 




l Qt 

2l-(Q*/Q ) 2 , 
Q t Qt 



Qt 



■ AW 



t- 



Qt 



e+lQi 



At 



(9) 
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where At is the constant time step and the components of increments AWt = 
W i+ At — W t , which are all independent real-valued random variables with 
mean and variance At. 



To solve equation (9) numerically, we need to give the values of Weissenberg 
number, Wi, and the maximum spring extension, Qn. Fo r convenience of com- 



parison with the experimental results of lSmith et al.l (119991 ). Weissenberg num- 
ber, decided by shear rate of the solvent and relaxation time of DNA molecules, 
is chosen in the range between to 80. Further, the contour length of A-phage 
DNA is chosen as the extensibility parameter. So, the maximum spring ex- 
tension can be expressed as Qq = (N — l)a, where N and a are the number 
of beads and the length of a rod of the corresponding bead-rod chain, respec- 
tively. All the parameters used in our computation can be found in Table 1. 
In the computation, the initial values of extension are chosen randomly as 
Gaussian distribution. 



Table 1. The value of parameters 



Parameters 


Values 


Parameters 


Values 


Weissengerg number Wi 


~ 80 


No. bead-rod chain N 


150 


Contour length of A-DNA L 


22yum 


Length of rod a 


0.147/zm 


Time step At 


10~ 5 ~ 10" 3 


Internal viscosity e 


0.001 



For the FENE dumbbells, there is a certain probability that the maximum 
allowed spring extension Q is exceeded for any finite time steps. To avoid this 
problem, we adopt the rejection method that all moves larger than a fixed big 
value are rejected. Because the extension that is very close t o the maximum 
allowed value will bring bad results, as proposed by lOttinger I (119961 ). we reject 
all moves which lead to a value of Q 2 larger than 



(l - y/At) Ql 



(10) 



and note that all the variables we use are in dimensionless form. 
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t is big enough? 



Output: Q 



Figure 3. Flow chart of the simulation 



After giving the initial position, the connector vector will be calculated by 
the algorithm stated above. Figure 3 shows the flow chart of the procedure. 
To obtain a good valuation of the mean extension, 1000 molecules or more 
should be considered, that is, the procedure will be performed repeatedly to 
yield more accurate mean value of the molecule extension. 



4 Numerical Experiments 



In this section, the computational method described in Section 3 will be ap- 
plied to the solution of equation (9) to obtain the extension distributions of the 
dumbbell of A-phage DNA moleculesfor different Weissenberg nu mbers. Sim - 



ulatio n results will be compared with the experimen tal data by ISmith et al 



(119991 ) and the simulation results of FENE model flHur fc Shaqfeh I . |2000|) 



as well as with the wor mlike spring model with hydrodynamic interactions 



(IJendreiack et al. I . l2002h . 



First, we study the extension's fluctuations along time. After the numerical 
experiments with the parameters listed in Table 1 for 2 x 10 6 dimensionless 
time steps with At = 10~ 4 , we obtained the time history of one DNA molecule 
extension at Wi = 6.3, as demonstrated in Figure 4. The fractional extension 
means the ratio between extension of one molecule with its contour length, 
Q/Qo- It can be found out from the raw data that there are large fluctuations 
in the extension of the molecules even at low shear rates. 
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Figure 4. Dimensionless extension versus time 



Next, we study the mean extension of 1000 molecules for different Weissenberg 
numbers. In the simulation, the molecular extension is projected in the x-y 
plane to match Smith's experimental data obtained by fluorescence microscopy 
with intensified video camera perpendicular to the flow-vorticity plane. Every 
simulation with different Wi has given a specific mean extension value after 
10 6 time steps when the average molecular extension appears steady. 



In Figure 5, we compare our simulation obtained with the FENE dumb- 
bell model accoun ting for internal viscosity with the experimental data b 



Smith et al.l (119991 ) . The dumbbell model without IV flHur fc Shaafeh [ 200 



and t he wormlike model with hydrodynamic interactions (IJendrejack et al. 



20021 ) are also superimposed in the plot. The classical FENE dumbbell model 



can not predict the extension result when Weissenburg number is high. Our 
dumbbell model and wormlike model have given better values at high Weis- 
senberg numbers. 
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Figure 5. Dimensionless extension of DNA molecules versus Wi 



When Wi = 0, that is, the shear rate of the solvent vanishes, the fluctuations 
of molecular extension are caused by the Brownian force bombarding on DNA 
molecules of small solvent molecules. When Wj increases , the e xtension reaches 



an asy mptotic plateau as found by ISmith et al.l (119991 ) and iHur &: Shaqfeh 



(120001 ). Our simulation results have better agreement with the experimental 
data compared to the other two simulations presented in Figure 5, especially 
when the value of Weissenberg number is high. We remark that this agreement 
with available experimental data is obtained with a relatively simple mathe- 
matical model where the internal viscosity effects were taken into account. 



Finally, we study the probability distribution of extension with the numer- 
ical methodology described in Section 3. As shown in Figure 4, from every 
snapshot of one DNA molecule we can get different values of extension. Fur- 
ther, every molecule at the same time shows a us different extension. Now, 
we investigate the probability distribution of the extension of DNA molecules 
for some specific Weissenberg numbers. After sampling 1000 DNA molecules 
during steady state, we plotted the probability distribution function (PDF) 
of the dimensionless extension in Figure 6. The PDF looks Gaussian in shape 
at small Weissenberg number (Wi = 1.3). Note further that when the Weis- 
senberg number increases, the peak of the PDF decreases and the distribution 
broadens. At large Weissenberg numbers (Wi = 76), th e PDF appears like 
white noise. For comparison, the experimental results by ISmith et al.l (119991 ) 
were superimposed in Figure 6. As seen, the simulation results have good 
qualitative and quantitative agreement with the experimental results. 
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Figure 6. The probability distribution function (PDF) of the dimensionless 

extension 

The demonstrated agreement with the experimental results is remarkable, 
given relative simplicity of the dumbbell model with IV and its computational 
efficiency. The model and the developed computational implementation pro- 
vide a tool for the analysis of extension and rheology properties. 



5 Conclusion 



In this paper, we focussed our analysis on the FENE dumbbell model with 
internal viscosity for the single DNA molecule in shear flow. The Brownian 
dynamics simulation was used to study the extension between the two beads 
of the dumbbell. The results are in agreement with the experimental data, 
which is especially important for high Weissenberg numbers. The probability 
distribution functions of the extension were also presented and the qualitative 
and quantitative agreement with available experiment data was demonstrated. 
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